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Abstract 

We examine the possibility of a confinement-deconfinement phase transition 
at finite temperature in both parity invariant and topologically massive three- 
dimensional quantum electrodynamics. We review an argument showing that 
the Abelian version of the Polyakov loop operator is an order parameter for 
confinement, even in the presence of dynamical electrons. We show that, in 
the parity invariant case, where the tree-level Coulomb potential is logarithmic, 
there is a Berezinskii-Kosterlitz-Thouless transition at a critical temperature 
(Tc = c^/Stt + (9(e'*/m), when the ratio of the electromagnetic coupling and 
the temperature to the electron mass is small). Above Tc the electric charge is 
not confined and the system is in a Debye plasma phase, whereas below Tc the 
electric charges are confined by a logarithmic Coulomb potential, qualitatively 
described by the tree-level interaction. When there is a topological mass, 
no matter how small, in a strict sense the theory is not confining at any 
temperature; the model exhibits a screening phase, analogous to that found 
in the Schwinger model and two-dimensional QCD with massless adjoint 
matter. However, if the topological mass is much smaller than the other 
dimensional parameters, there is a temperature for which the range of the 
Coulomb interaction changes from the inverse topological mass to the inverse 
electron mass. We speculate that this is a vestige of the BKT transition of 
the parity-invariant system, separating regions with screening and deconfining 
behavior. 
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1 Introduction 



One of the most intriguing features of a gauge theory is the possibihty of 
confinement. In a confining system, there are no "in" or "out" fields appearing 
in the asymptotic states which have color charges 0- As a result, all 
asymptotic states are singlets under the symmetry transformations in the 
color group. This generally occurs in one of two ways. First, as it is widely 
believed to be the case in four-dimensional quantum chromodynamics (QCD), 
the charged fields appearing in the bare action - quarks and gluons - are 
permanently confined into color singlet bound states - mesons and baryons 
- which make up the entire spectrum of asymptotic states. Second, as well aso 
the appearance of color singlet bound states, it is possible that the charges 
of bare fields are completely screened, so that they interpolate physical fields 
which occur in the spectrum but create color singlet states. This possibility 
has been raised for the electron field in the Schwinger model 0] and in two 
dimensional adjoint QCD where the bare quark mass is zero 0]. This situation 
is often referred to as "screening" rather than confinement. 

At finite temperature, the difference between a confined and deconfined 
phase is less evident than at zero temperature. There is no concept of 
asymptotic states and the quantitative observable features are thermodynamic 
variables and correlation functions governing the propagation of external 
influences. The commonly used test for confinement in a gauge theory at 
finite temperature is its ability to screen static external charges. The operator 
which probes electric screening is the Polyakov loop operator |^ ^ which is 
the trace of the path ordered exponential of the gauge field on a path which 
links the periodic Euclidean time Q 

P(f) =trP e'/o^^'^"^°("'"^) (1) 

The expectation value of this operator is the exponential of the free energy, 
F{x), which is required to immerse a classical fundamental representation 
quark source in the gauge field medium at the point x, 

(P(f)) = e-^(^^)/^ (2) 

If the expectation value is zero, corresponding to infinite free energy, this is 
interpreted as a signal of confinement. In this case, the system is not capable of 
screening the electric flux which is necessarily created with the external quark 
source and the electric flux takes up a configuration which has infinite energy. 
It is important that the infinite energy arises from an infrared, rather than 
ultraviolet divergence since the latter occurs even in non-confining theories 
and could be cured by introducing a fundamental cutoff. 

If the expectation value of the Polyakov loop is non-zero, and the free 
energy finite, this is interpreted as the system being in a deconfined phase. 
The electric fiux associated with the source is screened by the medium. 

In compact gauge theories, the Polyakov loop operator can be used as an 
order parameter for confinement in either pure Yang-Mills gluo-dynamics or 

^We use units where Planck's constant, the speed of hght and Bohzmann's constant are 
one. For a discussion of the path integral formulation of finite temperature gauge theory see 
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in a matter-coupled gauge theory when the matter is in either the adjoint rep- 
resentation, or some other representation whose degrees of freedom transform 
trivially under the center of the gauge group. In these cases, there is an invari- 
ance of the finite temperature path integral under gauge transformations which 
twist by an element of the center of the gauge group in the periodic Euclidean 
time. The operators in the action are all invariant and remain periodic (or 
anti-periodic in the case of fermions) under such a twisted gauge transform, 
but the operator P{x) is transformed to ZP{x) where Z is the central element. 
Thus, the finite temperature theory has an effective global symmetry whose 
transformations are the cosets of the set of all gauge transformations under 
which the path integral is symmetric modulo those which are strictly periodic 
in Euclidean time. Whether this symmetry is broken or not is a well defined 
question. Spontaneous breaking of the symmetry is related to deconfinement 
p|, P] and P{x) is an order parameter. This order parameter has been par- 



ticularly useful in characterizing the nature |p|-|ll] and quantitative aspects 



T2|-|n| of the confinement-deconfinement phase transition in a wide array of 



pure gauge theories. 

When the gauge theory is coupled to matter fields which transform non- 
trivially under the center of the group, the symmetry is broken explicitly and 
its realization can no longer be used as a probe for confinement. An example 
is QCD with quarks in the fundamental representation of SU(3). In that case, 
the question of distinguishing a confining and non-confining phase of the finite 
temperature gauge theory is more sophisticated |]T5 |. 



Recently it has been argued |T6|, that the Abelian analog of the Polyakov 
loop operator, 

P,{x) = e'^^o^^'^^ (3) 

could be used as an order parameter for confinement in Abelian gauge theories, 
even in the presence of dynamical charged particles. The requirements are that 
the dynamical charged fields in the gauge theory must have charges which are 
integral multiples of some basic charge, which we denote by e. Also, a technical 
requirement is that the charged fields have a mass gap and that the field theory 
has a finite ultraviolet cutoff. In the limit where the cutoff is removed, it is 
usually necessary to define the loop operator by multiplicative renormalization. 
In order that (^) be an order parameter, the charge e appearing there must 
not be an integer multiple of the basic charge: e ^ e ■ integer. 

The expectation value of the operator (j^) measures the response of an 
electrodynamic system to placing a classical incommensurate charge e at point 
X. The quantity 

F,{x) = -T In {P,{x)) (4) 

is the free energy of the system in the presence of the classical charge (minus 
the free energy when the charge is absent). If this free energy is finite, the 
system is not confining. If the free energy is infinite, this implies that the 
expectation value (|^) must be zero. This means that it takes an infinite amount 
of energy to immerse a classical charge in the system, implying that it is in a 
confining phase. If the charge in the loop operator were a commensurate, rather 
than incommensurate one, its electric field could be screened by producing 
a finite number of dynamical charged particle-antiparticle pairs, using the 
appropriate number of particles or antiparticles to screen the external charge 



3 



and allowing the remaining dynamical particles to escape to infinity. This 
process would take a finite amount of energy. For this reason, we expect that 
the loop operator with a commensurate charge would always have a non-zero 
expectation value. At zero temperature, the pair production would take a 
threshold energy of the mass of the particles produced. At finite temperature 
the thermally activated particles are already present so that screening of this 
sort takes a small amount of energy. 

Further information can be obtained from the correlators, 

Fe,..e„(fi, . . . = -Tin P,M)^ (5) 

which give the electrostatic energy of an array of charges Cj situated at points 
Xj, respectively. For example, the two-point correlator gives the effective 
potential between a positive and negative charge, 

V,,_,{x,y) = Fe,_e(f,y) = -T\n{P,{x)P^M) (6) 

If the clustering property of this expectation value holds, vanishing or non- 
vanishing of the expectation value of a single loop operator is related to the 
asymptotic behaviour of the potential: if 

lim Ve-e{x,y) = oo (7) 
the two point function has the behavior 

^ lim (Pe(x)P_e(f)) = (8) 

and the cluster decomposition property implies that the expectation value of 
the single operator should vanish 

(Pg(f)) = (9) 

This characterizes confinement. If, on the other hand 

lim Ve-e{x,y) = finite constant (10) 

\x~y\^oo ' 

the potential is not confining, the two-point correlator has the behavior 

lim {Pe{x)P-e{y)) = constant , (11) 

\x-y\^oo 

clustering implies that 

{Pe{x)) = constant (12) 

and an isolated external incommensurate charge has finite energy. 

The expectation value of the loop operator in (|) is governed by a particular 
discrete global symmetry, Z, isomorphic to the additive group of the integers, 
which appears in the euclidean path integral at finite temperature. To see its 
origin, consider the euclidean finite temperature functional integral expression 
for the partition function 0, 

Z[T]= I n [dA^{T,x)dij{T,x)di!{T,x)]e-^^^^'^^^ (13) 
x,T6[o,i/r] 



4 



where the action is 

S[A, tp, ^] = j^J^ dr J dx (^^F^A^, x) + tP{T, x) {-f^D^ + m) ^(r, x)^ (14) 

with = — ieA^ and F^^, = V^A^ — V^A^. The boundary conditions in 
time are periodic for the photon, 

A^{t = 1/T,x)=A^{0,x) (15) 

and antiperiodic for the electron, 

ij{T = 1/T,x) = -4j{t = 0,x) , V^(r = 1/T,f) = -V5(r = 0,x) (16) 



The path integral (|T^) is symmetric under gauge transformations, 

A;(r, x) = A^{t, x) + V^x(t, x) (17) 

^'(r,x) = e^"'^("'^V(r,^) , ^'{t,x) =iIj{t,x) e-'''''-^''^^ (18) 
when the gauge transformation function x(r, x) has periodic derivatives, 

V^x(r = 1/T, x) = V^xir = 0, x) (19) 

and it is periodic up to an integer multiple of 27r/e, 

x{t = 1/T, x) = x{r = 0,x) + 27m/e , neZ (20) 

The group of all gauge transformations modulo those which are strictly periodic 
is Z, the additive group of the integers. This is a global symmetry. The 
Polyakov loop operator transforms non-trivially under the coset when its 
charge is not an integer multiple of the electron charge, 

P^{x) = Pg(f) . e^™^/^ (21) 

It can therefore be used as an operator to explore the realization of Z in 
the statistical model specified by the path integral (0). If the symmetry is 
unbroken, the loop operator averages to zero and the system is in the confining 
phase. If it is spontaneously broken, the loop operator can have a non-zero 
expectation value. The system is then in a non-confining phase. We shall 
discuss two such phases. One is the Debye plasma phase which is characterized 
by the exponential decay in the asymptotic behaviour of the electric two- 
particle potential, related to the Debye screening length of the plasma. It 
is also related to the mass of elementary excitations in a lower dimensional 



theory [^]. The other is what is termed a "screening phase". This phase 
also has an electic screening length, as well as magnetic screening. The main 
difference between these two phases is in the temperature dependence of the 
screening length. In the screening phase, the screening length persists at zero 
temperature, whereas the Debye mass vanishes at zero temperature. 

The Z symmetry has a physical interpretation in terms of the charge of 



physical states |]19[. In path integral quantization, as we shall discuss in Sec- 



tion III, the temporal component of the gauge field Aq arises as a Lagrange 
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multiplier to enforce gauge invariance. The projection operator which guar- 
antees gauge invariance in the construction of the path integral is obtained 
by exponentiating the generator of infinitesimal gauge transformations and 
integrating over all gauge transformations, 

-p = _1_ J [dAoix)] / '^"^ {VAo-E-Aoe:^H:) (22) 

and the density matrix for the Gibbs ensemble is 

Q-H/T 

P = P^ (23) 

The Z transformation 

^o(^) ^o(^) + 27mT/e (24) 

results in 
where 

Q = e[dx : ^\x)^/j{x) : (26) 



is the (normal ordered) electric charge operator. If the Z symmetry is not 
broken, all physical states with non-zero weight in the thermal ensemble have 
quantized charges, 

-H/T -H/T 

^ ' Z[T] ~' Z[T] 

If the Z symmetry is broken, there are states contributing to the thermal 
ensemble which have non-quantized charges. Intuitively, these can occur in a 
deconfined theory since the long-ranged electric fields accompanying arbitrarily 
diffuse charge distributions would have finite energy. Such states would not 
be allowed in the confining phase. Note that, in this operator picture, the 
Z symmetry is not a symmetry transform of the density matrix in the usual 
sense. In fact the existence of this symmetry is related only to the question of 
whether exp(27ri Q/e) is the unit operator when operating on the states that 
have non-zero weight in the density matrix. 

At T = 0, and for the physical value of the electromagnetic coupling 
constant, 3+1-dimensional electrodynamics does not exhibit a confining phase. 
It is in the deconfined Coulomb phase at zero temperature and forms a Debye 
plasma at finite temperature and density. There is a conjecture that if the 
electric charge of QED could be increased to some critical value, there would 



be a phase transition to a chiral symmetry breaking and confining phase |20]. 
It is reasonable to expect that the resulting strongly coupled system would 
have a confinement- deconfinement transition at some finite temperature. 

QED in 1+1 dimensions with a massive electron, i.e. the massive Schwinger 
model, is confining and the Z symmetry is not broken. In the case of the mass- 



less Schwinger model, the Z symmetry is spontaneously broken [21, 19]. This 
spontaneous breaking of Z is interpreted as screening, rather than deconfine- 
ment. A similar situation appears in 1+1-dimensional QCD with massless 



quarks P^l ffl. With massive quarks, the question of confinement and the 
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confinement-deconfinement transition in the large limit of 1+1-dimensional 
QCD both at finite temperature [ P5| , ^ and at zero temperature |^ has 



been examined by several authors. 

When the electron has a mass, parity invariant quantum electrodynamics 
(QED) in 2+1 dimensions is believed to be a confining theory. The tree level 
Coulomb interaction varies logarithmically with distance. Its entire spectrum 
is bound states, although the bound states can have arbitrarily large sizes. 
The perturbative self-energy of the electron has a logarithmic infrared infinity 
P7| . Also, when the number of electron fiavors is small enough, confinement, 
accompanied by chiral symmetry breaking is believed to persist in the limit 
as the bare electron mass is put to zero [^. In a previous paper by three 
of us it was argued that 2+1-dimensional QED has a phase transition from 
the confined to a deconfined phase at a critical temperature [|l^. The phase 



transition is of Berezinsky-Kosterlitz-Thouless [^9|, ^ (BKT) type with non- 
universal, temperature dependent critical exponents. This is the standard 
phase transition of the Coulomb gas, which is a reasonable characterization of 
the thermal state of QED in the limit where the density of thermally excited 
particles is low. This is the case when the mass of the particles is greater than 
the temperature. The Polyakov loop operators have power-law correlators 
in the confining phase. In the deconfined phase there is an electric, Debye 
mass which makes the Coulomb interaction short-ranged. There is a universal 
quantity associated with the BKT phase transition, the bulk modulus of spin 
waves in the massless phase. This predicts the power law behavior of Polyakov 
loop correlators in the confined theory. 

It is interesting to ask what happens in the case where the three dimensional 
gauge theory is not parity invariant, but the action contains a topological mass 
term for the photon. Of course, one would expect that the resulting photon 
mass 



|32| provides an infrared cutoff and thereby removes the long-ranged 
interactions that are responsible for confinement. However, if the topological 
mass is very small, so that the confinement scale is much larger than the photon 
mass, we expect that some effects of confinement persist at short distances. 
To understand what effect topological mass has on the Z symmetry. 



consider the Chern-Simons (CS) action in Euclidean space 

Scs = ^\ J^^^ dr J dxeoij {2AoViAj - AAj) (28) 
The gauge transformation 

Ao(r, x) ^ Ao(r, + ^ (^^^'") ^^9) 

where e is the basic unit of charge of all matter fields, results in the change of 
the action 

Scs Scs + ineg (30) 
where g is the magnetic charge, 

g = J dxeoijViAj (31) 

When, as it happens on an open space such as TZ^ which we consider here, the 
magnetic charge is not quantized, the Chern-Simons action is not invariant 
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under the gauge transformation unless n = 0. The presence of a Chern-Simons 
term in the action then would break the Z symmetry explicitly. 

On the other hand, if the space is compact, or if we impose boundary 
conditions on an open space, such as the plane, so that the gauge field can be 
stereographically projected onto a compact space, the magnetic charge should 
obey the Dirac quantization condition, 

eg = 27rk (32) 

where k is an integer and e is the basic unit of charge of the matter fields. 
In this case, the Chern-Simons term is invariant under gauge transformations 
with non-zero winding number. The global symmetry is Z, even in the absence 
of matter. 

The Polyakov loop operator for a basic charge (remember that, because of 
the presence of monopoles, the charge is quantized) transforms as 

gie/g^^drAo(r,x) ^ ^je^/K _ ^ie J^^^ dTAo{T,x) j-ggj 

We denote the coefficient of the Chern-Simons term as 

so that the Polyakov loop operator transforms by the phase e^'^**/^. Depending 
of the ranges of p and q, there are three possibilities: 

• li p/q is an irrational number, then the symmetry is Z. In finite 
volume this will imply that all charges are confined and only neutral 
configurations are allowed. 

• If p and q are integers so that p/q is rational then the symmetry is 
the finite cyclic group Zp. In this case, charges which are not integral 
multiples of p are confined. 

• If g is an integer and p — 1, there is no symmetry. 

The latter condition is compatible with Gauss' law which, as we shall see 
in Section III, relates the total charge and magnetic flux of a quantum state 
on a compact space as 

2 

+ (35) 

Since g — ^ ■ integer and Q = integer • e, it is necessary that p and q are 
integers. Furthermore, the basic state in this system has electric charge p ■ e 
and magnetic charge q - —. 

In fact, we could think of the effective symmetry of the Euclidean path 
integral as just enforcing the global constraints on the charges which are 
contained in Gauss' law. We shall examine the question of whether this 
symmetry survives in the infinite volume limit in Section III. 

In order to analyze symmetry breaking, the properties of correlators and 
other dynamical questions in parity invariant 2-|-l-dimensional QED, we shall 
compute the effective action for the Polyakov loop operator in Section II. This 
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method was advocated in the seminal work Svetitsky and Yaffe in the context 
of lattice gauge theories [0. One first fixes the static temporal gauge, 



^Ao(r,f)=0 (36) 

and integrates all the degrees of freedom of the gauge theory except for Aq 
which is associated with the Polyakov loop. This generates an effective theory 
for the order parameter which, for an initial d+ 1-dimensional gauge theory is a 
d-dimensional scalar field theory with variable Af){x). This theory exhibits the 
global Z symmetry explicitly. Since, at finite temperature, fermions always 
have a mass gap, integrating out all the other degrees of freedom in the 
original theory produces only short-ranged interactions in the effective theory. 
Consequently, the critical behavior oi d + 1-dimensional finite temperature 
QED is that of the d-dimensional effective local field theory. By studying the 
effective action, we are able to characterize the type of the phase transition 
and discuss the associated critical behavior. 

In Section III we shall consider 2+ 1-dimensional pure QED with a Chern- 
Simons term in the presence of an array of external charges. We consider 
the case where the space is compact and discuss some boundary problems 
related to the presence of the Chern-Simons term on a compact surface. We 
discuss canonical quantization and the construction of the functional integral 
representation of the partition function Z[T]. We then find the exact effective 
action for Ao by integrating all the spatial components of the gauge field. We 
use this effective action to study the realization of Z symmetry in the infinite 
volume limit. 

In Section IV, we use a variational method to examine the behavior of both 
parity invariant and topologically massive QED in infinite volume. We confirm 
the existence of the BKT phase transition in the parity invariant theory. We 
also find indications of a similar transition in topologically massive QED. We 
speculate that the latter transition, although not strictly speaking a phase 
transition, separates two regions of the system with distinct physical behavior: 
a low temperature screening phase where there are charged particles bound into 
neutral bound states, as well as free neutral particles, and a high temperature 
deconfined phase where the bound states are absent and the charges of the 
particles are Debye-screened. Section V is devoted to a discussion of our 
results. 

In this paper, we do not address the interesting and controversial question 
of whether domain walls exist between regions with different "orientations" 
of broken Z symmetry. The existence of these in non-Abelian gauge theories 
has been a subject of much discussion ^ 0, |3^. If they did exist in 



QED, they would be very interesting and perhaps observable objects. Detailed 
analysis of this possibility is still an important problem. 

It is interesting that in this paper we find a phase transition which is 
accessible to perturbation theory. This is a property of the critical line of 
BKT transitions; one end of the line is in a perturbative regime. To our 
knowledge, this is the only situation where one can study a confinement- 
deconfinement transition without the aid of numerical simulations. On the 
other hand, confirming the existence and properties of the transition which we 
discuss by numerical simulations would be a most worthwhile project. 
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2 Effective action in parity invariant QED 



We shall consider QED in 2+1 dimensions. As it is well known, in 2+1 
dimensions, the minimal, two-component Dirac fermions violate parity [^, 
0, so that if included in the action, they can generate a parity violating 
topological mass for the photon by radiative corrections ^ . In this section 



we shall study the case where the electron has mass but the photon is massless. 
For this purpose, we shall use parity invariant four- component fermions which 
are obtained by dimensional reduction of the 3+1- dimensional Dirac operator. 
The resulting model has two species of massive two-component fermions where 
the mass terms have opposite signs. The Euclidean action is 



S= (fx 



^F2^ + ^(7- (V-zeA) +m)^ 



(37) 



We shall find the effective action for the Polyakov loop operator (^. At finite 
temperature, it is possible to use a gauge transformation to set the temporal 
component of the gauge field, y4o(r, x), independent of the Euclidean time, r. 
In that gauge, 

p-{x) = e'-'^«("^)/^ (38) 

and the effective field theory for this operator is the effective action for the 
static field Aq{x). This two-dimensional effective field theory is obtained by 
integrating the other degrees of freedom from the path integral 

^-SM] ^ f [dAdtpd^] (39) 



(The functional integral defines the effective action only modulo a temperature 
and volume dependent but Aq independent constant.) The integral over the 
Grassman-valued Fermi fields yields the determinant of the Dirac operator, 
and 

— * — » 

S^f,[Ao] = J dx^^^^;^^ - In J[dA] det[7-(V-zeA) + m] e'/o^^'^"''"' ^^V^ 

^ (40) 
The first term in (|40|) is the tree level action and the second term contains 
all quantum corrections. The remaining functional integral requires additional 
gauge fixing. The Z symmetry is a periodicity of the effective action under 
the field translation 

S,s[Ao] ^S,s[Ao + 27rnT/e] (41) 

The effective action, S'efT[Ao], is non-local and non-polynomial. It can only 
be regarded as a local field theory when the momenta of interest are much 
smaller than the mass gaps of the fields which have been eliminated. In that 
case the effective action has a local expansion in powers of derivatives divided 
by masses. 



At finite temperature the action (|37|) contains three parameters with 
the dimension of mass, the electron mass m, the gauge coupling e^, and 
temperature T. The loop expansion is super-renormalizable [^1], In fact, 
with a gauge invariant and Euclidean Lorentz invariant regularization, it has 
no divergences whatsoever. This means that the effective theory that we obtain 
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does not contain an ultraviolet cutoff, and its only dimensional parameters are 
e^, m and T. 

The dimensionless parameter which governs the accuracy of the loop 
expansion is the smaller of e^/m and e^/T. Also, in order for the local, 
derivative expansion of the effective action to be valid, it is necessary that 
either the electron mass or the temperature be larger than the momentum 
scales of interest. Then, the larger of m or T acts as the ultraviolet cutoff of 
the effective field theory. 

At the tree-level, where the effective action is approximated by the first 
term on the right-hand-side of (^0]) only, we can easily compute the correlator 
of Polyakov loop operators, 

ieAoix) ^-^eAo{0)\ _ coUSt. J] 1^1"''^^"^ (42) 



e ' e 



/ tree 



It has a scale invariant form with a temperature dependent exponent reminis- 
cent of the correlators in the Gaussian spin- wave theory This is a 



result of the marginally confining nature of the logarithmic Coulomb interac- 
tion. The behavior is between that of a confining theory where the correlator 
would exhibit the clustering property and decay exponentially at large dis- 
tances and the deconfined theory where it would approach a constant. It is 
interesting to ask how quantum fluctuations would modify this result. We shall 
argue in the following that, at low temperatures the behavior ( ^2] ) is qualita- 
tively, though not quantitatively correct. At high temperatures, the correlator 
approaches a constant at large distances and the Z symmetry is broken. 

We shall compute the effective action for Aq{x) in the one- loop approxima- 
tion and in an expansion in powers of derivatives of Ao{x). The leading terms 
are 

S,s[Ao] = J d?x {Z{m, eA,/T)^VA^ ■ VAq - V{m, e^o/T)) . (43) 

Here V is the effective potential for Aq. Z — 1 arises from the expansion of 
the temporal components of the vacuum polarization function to linear order 
in -V^. 

noo(^ = 0, k^) = noo(0, 0) + fc2(Z - 1) + . . . (44) 

To one-loop order, the effective potential, V{m,eAQ/T), is obtained from the 
fermion determinant in the constant background Aq, 

V{m, eAo/T) = — ^ log det((-^9o - eAof - + m^) (45) 
(Vol.) 

The fermions have anti-periodic boundary conditions in the compact time. To 



regularize the determinant we consider the ratio ^ 



/I m det((-zgo-eAo)^-V^ + m^) 

= det(-ag-V^ + m^) ^^'^ 

The antiperiodic boundary conditions lead to the expression 

n.k 
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The product over integers can be evaluated explicitly as |^ 

sin2(eAo/2r) 



A(m,eAo/r) = n 1 

k - 



cosh2(Afc/2r)_ 



= Y[A^^{m,eAo/T) 

k 



(4J 



where = fc^ + are the eigenvalues of the operator — + m^. In the 
infinite volume limit the 11^ in Al'^^eAo/T) gives rise to an integral on k 
in log A. The result (^) holds in any dimensions, but the integral on k 
can be performed analytically only in two dimensions. In one and three 
dimensions this integral can only be done for m = in which case it gives 
simple polynomial expressions. In the limit m = 0, the effective potentials for 
^0 have been discussed in |16|. In two dimensions we obtain 

'•+°° d'^k 



V{m,eAo/T) 
T2 



f27rl 



logA,^(m,eAo/T) 



TT 



m 
T 



Li2{ e-'"/^, eAo/T + vr) + Li^{ e"™/^, cAq/T + 



where 



and 



rr 

Li2{r,6) = - dx\n{l -2xcos6 + x'^)/2x 
Jo 



dxLi2{x, 0)/x 



(49) 

(50) 
(51) 



are the real parts of the dilogarithm and trilogarithm according to the 
convention of Ref. [^. As one can see from the definition of Li2{r,6) and 
Li3{r,9), the effective potential is explicitly invariant, as expected, under the 
Z symmetry, cAq/T cAq/T + 27m. 

Computing the temporal components of the vacuum polarization function, 
noo(p, ^o)? ill an external constant Aq field and keeping only the term which 
contributes the leading order in derivatives to the effective action, one obtains 



Z{m,eAo/T) = 1 + 



127rm 



1 — m 



d 



sinh m/T 



dm J cosh m /T + cose Aq/T 



(52) 



We will find it convenient to use the harmonic expansion of the effective 
potential (^), 



V{m,eAo/T) 



(53) 



There are two limits in which we can obtain more analytic information 
about the effective potential. In the regime m » T, T/m and e^/m are small 
and e^/T is unrestricted. The higher harmonics in the effective potential are 
exponentially small perturbations to the leading term. 



Vim,eAo/T) 



Tm 



TT 



-m/T 



cos(ey4o/T) 



(54) 



which is the sine-Gordon potential. The effective field theory in this limit is 
thus the sine-Gordon model, 



S. 



{m»T) 
eff 



[A, 



dx 



Tm 



—VAo ■ VAo - 

Zl 71 



-m/T 



cos(ey4o/T) 



(55) 
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This model has a phase transition corresponding to the BKT |30| transition 
in a two-dimensional classical Coulomb gas. The critical behavior associated 
with this transition has been studied extensively 



44, 45, 46, 47, 4g, 49 



Wiegmann |^ and Amit et al. |^9[ showed that in the sine-Gordon model any 
perturbation of the type cos(?t,/50) to the sine-Gordon potential acos(/50)//3^ 
is irrelevant to the critical behavior of the model. They showed that the 
scale dimension of cos(n/30) is 2n^ at the critical point. Consequently for 
n > 1 these operators have scale dimension greater than two and they are 
indeed irrelevant. Thus we conclude that the critical behaviour in the regime 
m >> T, is identical to that of the two dimensional sine-Gordon potential 
of eq.O). 

The BKT transition occurs at a line of critical points. In the sine-Gordon 
model with potential 

acos(/50)//3^ 
the critical line begins at the point 



(56) 



This critical point was originally found by Coleman in his discussion of 
bosonization and the correspondence of sine-Gordon theory with the massive 
Thirring model The line of critical points in the sine-Gordon theory 

corresponds to a line of critical points for the confinement-deconfinement 
transition in QED which can be drawn for example in the e~™'/^,e^/Tj 
plane. Comparing eq. ( ^4]) with the sine-Gordon potential we see that the 
QED critical line starts at 



m 



eVT =(0,87r) 



The critical temperature corresponding to this point is 



(57) 



(5^ 



The renormalization group was used to study this phase transition, originally 
pn| and Wiegmann and later improved to higher order by 
The flow diagram is depicted in Fig. 1. There are 



by Kosterlitz 
Amit and collaborators [H9 



three regions: The high temperature, deconfined region III, where the model 
is asymptotically free. This is the case which was analyzed by Coleman ||44|| . 



The low temperature region I is confining and has a line of infrared stable fixed 
points at m = cxD, corresponding to c = 1 conformal field theory. In Region II, 
the model is deconfined and is neither asymptotically nor infrared free. The 
separatrix between regions I and II is the line of BKT phase transitions. 

To compute the leading correction to Tcrit. due to a finite (but still large) 
value of the fermion mass, one has first to renormalize the field Aq, according 
to 

Al'"" = AoZ{m,OY/^ (59) 

This renormalization changes the sine-Gordon parameter (3 in the argument 
of the cosine. The one-loop calculation of the effective action results in the 
correction 



(m»T) 



^7r(l + eVl27rm + . . .) 



(60) 
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As we have shown, the BKT phase transition in QED is a confinement- 
deconfinement transition. In the spin wave plus Coulomb gas description of the 
XY-model the BKT phase transition corresponds to a binding- unbinding 
of vortices. In QED it has the obvious analog of a binding-unbinding transition 
for charged particle-antiparticle pairs. In the deconfined phase, Aq fluctuates 
near one of the minima of the effective potential, 

{Ao) = 2TxnT/e (61) 

In a semiclassical analysis, this expectation value contributes an imaginary 
chemical potential for the electron. However, this chemical potential can 
be absorbed by shifting the Matsubara frequency by n units. Thus, in a 
semiclassical analysis, the thermodynamics in the deconfined phase does not 
suffer from the difficulties of negative entropy and imaginary thermodynamic 
potential that affect the meta-stable Zj^ phases of QCD ||5y, |5T1] . 

The other limit where the one-loop result is simple is the high temperature 
limit, where T >> m, e^. In that limit we must be careful to study the degree 
of freedom which is periodic 0. For this, we define the field 

a{x) =eAo{x)/T (62) 

so that the effective action for the field a{x) has the periodicity 

a(x) a{x) + 27m (63) 

The effective action in the high T limit is 

sZ^^""\a]= I ^Va-Va+ ^Lz3(l, a + 7r) I (64) 

Large T is the semi-classical limit for this theory and a must fluctuate near a 
minimum of the effective potential. In this expected, the Z symmetry 

is spontaneously broken, corresponding to deconfinement. 



3 Maxwell-Chern-Simons theory on the sphere 

It is interesting to ask what happens in three-dimensional QED when it is not 
parity invariant. In this case, the gauge field can have a topological mass term 



3T| and, naively, one would expect that confinement is not an issue, it is simply 
absent. In this Section, we shall consider the properties of finite temperature 
Maxwell-Chern-Simons theory when the space is the 2-sphere. We shall find 
that there is an analog of the Z symmetry, which exists and has interesting 
properties even in the absence of matter fields. The symmetry enforces a kind 
of topological confinement which arises from Dirac's quantization condition 
for the magnetic field of the monopole. For completeness, we shall also give 
a careful treatment of the finite temperature path integral in this case. The 
Minkowski space action for Maxwell-Chern-Simons theory with Wilson-loop 



sources is 



S = Jd'x {-^F^.F^'^ + ^e^-M^9,A,| + E / (65) 



We disagree with the discussion on this point in ref. Ilf 
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The spacetime is a product of time and S"^. A similar construction can be 
carried out where the space is a product of and any Riemann surface with 
some additional complications ||52|| . 

Since the space is compact and we shall consider the situation where the 
total charge of the external charge distribution is not zero, the Gauss' law 
constraint will force us to consider the case where there is a non-zero magnetic 
charge 

' dA = g^O (66) 



/52 

In this case, the spatial components of the gauge field are not globally defined 
functions on the sphere but are rather components of a connection on the 
monopole line bundle. It is well known that, in this situation, extra topological 
terms are needed to make the Chern-Simons term well-defined |^, This 
arises from the fact that the density which one integrates to get the Chern- 
Simons term is not gauge invariant, but transforms by an exact form under 
gauge transformations. This makes the integral of the Chern-Simons three- 
form sensitive to the coordinatization of the monopole line bundle. 

This sensitivity can be seen by the following argument. We decompose the 
gauge field into its spatial and temporal part, 

Af.dx'^ = Aodt + Aidx' = A^dt + A (67) 

We construct the monopole line bundle by considering a set of coordinate 
patches, {-P/t}, which cover the 2-sphere, 

Ufc Pk = (68) 

and denoting the gauge field in the fc'th patch as A^ . The gauge field in 
different patches are related by gauge transformations, 

A^ -A^ = dx''" (69) 

where x''^ is a function which is defined on the intersection region P^UP^. The 
integral over S'^ is defined using a partition of overlapping patches, so that each 
point of 5*^ is integrated only once. Gauge invariant quantities such as dA or 
dAdA are not sensitive to details of the choice of coordinate patches. Likewise, 
the integral of gauge invariant quantities does not depend on the positions of 
the patches. However, the integral of a non-gauge invariant density, such as 
the (naive) Chern-Simons term 

AodA'' + A'^dAo - ^'"^^'"^ (70) 

depends on the position of the patches. For example, the contribution to 
Gauss' law density arising from the Chern-Simons term is obtained by taking 
a variational derivative of the Chern-Simons term by Aq. To do this, it is 
necessary to translate Aq —>■ Aq + SAq, to isolate 6Aq by integrating by parts 
and then to identify the functional derivative as the coefficient of 6Aq under 
the integral. With this procedure one obtains the charge density 

jQ = 2dA + J2 f 6{x - n Pp)dx'^ (71) 
kp •'P^^Pp 
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This has the unappeahng feature that some charge hves on the arbitrarily 
chosen transition regions (in fact, on the arbitrarily chosen boundaries between 
patches which are inside the transition regions). The surface terms which must 
be added to the Chern-Simons term cure this difficulty. In the present case, 
they are integrated on the intersection regions of coordinate patches and cancel 
the terms obtained when the second term in the above naive Chern-Simons 
term is integrated by parts to obtain the expression 




This version of the Chern-Simons term will be sufficient for our purposes in 
the following. 

3.1 Canonical quantization 

We shall examine the electrostatics of an array of classical charges. In this 
case, the source term in the action has the form 

J2e,JdtAo{x,,t) (73) 

i 

In particular, we are interested in the free energy of this system as a function 
of particle positions. We shall not require global neutrality of the charge 
distribution. However, the consistency of the monopole bundle will force us to 
use charges which are integer multiples of a basic charge, e, compatible with 
the Dirac quantization condition. 

If the total charge is non-zero, the gauge constraint, which is obtained by 
taking a functional derivative of the action ( p5| ) by Aq is 

V ■ E + kB + Y, eiSix -Xi) r^O (74) 

i 

The electric field is gauge invariant and must be a globally defined vector field 
on S'^. Therefore, the integral of the divergence of the electric field over the 
space 5*^ must vanish. The integral of Gauss' law then implies that, when the 
total electric charge is not zero, there is also a non-zero magnetic fiux, 

Kg + J2e, = (75) 

i 

It is convenient to separate the effect of the magnetic background field by 
decomposing the gauge field into a classical time-independent part containing 
the monopole field and a time dependent part with no overall magnetic fiux 
and which is allowed to have quantum fiuctuations, 

Ai{x, t) AmAx) + Ai{x, t) . (76) 

In eq.(|7|) 

/ VxAM = g (77) 
and Am is defined in such a way that the classical magnetic field is constant, 

Bm = V X Am = g/^TTR^ (78) 
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(with R the radius of S*^) so that 

Jb = o 

and 



J BmB = . 
Substituting ( [76D into the action yields 

S = Jdt J^ |i (i, - V.Aof - ^{Bl, + B^) + t,A,BM + f^A^B 

X i - 1^ (Am X A) I + X y dtCiAoit, Xi) (79) 



where, now all variables are single- valued, globally defined functions on S"^ and 
the monopole gauge field Bm is a classical variable. The only remaining multi- 
valued term / ^^(Am x A) is a total time derivative term and therefore is not 
important for the canonical quantization which we shall do in the following. 
Canonical quantization proceeds by identifying the canonical momenta 



Ho ~0 (80) 

— ( 

2 



Hi = Ai-ViAQ + -eijAj (81) 



The first relation is a primary constraint. The Hamiltonian is 

' 1 /„ K .V 1 „o 1 



and the canonical commutation relations are given by 

[A,(f),n,(f)] = i6,,6{x-y) (83) 
[Ao(f),no(^)] = i5{x-y) (84) 

The gauge constraint arises as a secondary constraint from requiring that the 
primary constraint Hq ~ is time-independent 

G{x) = V ■ n(f) + -B{x) + kBm + ^i^i^ -Xi) r^O (85) 

The operator G{x) generates time- independent gauge transformations and 
commutes with the Hamiltonian which is invariant under the gauge transfor- 
mation 

A{x) A{x) + |l(f), J dy x(y)G'(y)| = A{x) - Vx(x) 
n{x)^n{x) + S^U{x),Jdyx{y)G{y)^ = n(f) - |v*x(^) (86) 

where V* = ejjVj. 

The dynamical system with Hamiltonian (|82D and canonical commutator 
is internally consistent and can be quantized as it is (with the subtlety 
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that the ground state is not normahzable) . The primary constraint, ([80| ) is 
solved by imposing the auxihary gauge fixing condition 



Ao{x) ~ (87) 



and thereby ehminating both Aq and Ho. The Gauss' law constraint (p5| ) 
must then be imposed as a physical state condition. Since the operator G{x) 
commutes with the Hamiltonian, it can be diagonalized simultaneously with 
the Hamiltonian. Then the simultaneous eigenstates of the Hamiltonian and 
G{x) which are in the kernel of G{x) are chosen as the physical states, 

G{x) \ physical state >= . (88) 

We can form a projection operator onto physical states by considering the set 
of all gauge transforms. Formally, 

V = ^ j [dx{x)] exp [i J^^ x{x)G{x)^ . (89) 
This operator satisfies the property 

(This is a formal statement due to the infinite volume of the gauge group. This 
is the same difficulty which appears in the normalization of the states.) 

The form of Gauss' law indicates that the gauge symmetry is realized 
in a projective representation. For example, if we represent the canonical 
commutation relation in the functional Schrodinger picture, where states are 
wave-functionals "^[A] of the classical field A, and the canonical momentum is 
a functional derivative operator. 

Then a physical state which obeys the gauge condition 

G{x)^ptys[A] = (91) 

gauge transforms as 

^phys[^- Vx] = e*«/x(^M+^/2)+^E.->x(-^.)vi>p,yjl] (92) 

In the next subsection, we shall find the functional integral representation of 
the thermodynamic partition function. 

3.2 Functional integral representation of the partition 
function 

In this subsection, we shall discuss the derivation of the functional integral 
representation of the thermodynamic partition function. It is obtained by 
taking the trace over physical states of the Gibbs distribution operator 

p= e-^/^ (93) 
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where H is the Hamiltonian and T is the temperature. (We work in a system of 
units where the Boltzmann constant, the speed of hght and Planck's constant 
are equal to one.) We shall consider the unconstrained space of states which 
represent the canonical commutation relation (83) and insert into the trace a 
projection operator which projects onto the physical states: 

Z[T] = E < ^1 e-^'/^Pls > (94) 

s 

Explicitly, this can be written as 

Z[T] = — j [dx{x)] j [da,{x)] < a\ e"^/^ j x{BM+B/2)+iY,e,x{s^)^(^j^^^ > 

(95) 

where we have taken the trace using the eigenvectors of the "position" operator 

Ai{x)\a >= ai{x)\a > (96) 



The integration over x iii (ED projects the trace onto gauge invariant states. 
Since the Hamiltonian is gauge invariant, it is sufficient to insert the projection 
operator once. The field x{^) is proportional to the temporal component of 
the gauge field, Aq{x), in the time-independent gauge 

x{x) = Aq{x)/T . (97) 

We are particularly interested in deriving an effective action for the gauge 
function Ao{x). It is defined by 

^-sMAo] = ^ / [da] < d\ e-^/^ e'?/^»(^«+^/2)|^^ > (gg) 

Bai 

where we have omitted the external charges and we have summed over 
magnetic monopole number. The subsequent integration over Ao{x), which is 
needed in order to obtain the partition function of the system in the presence 
of external charges, will enforce Gauss' law. In particular it will project onto 
the sector with the correct magnetic charge. 

The partition function of the system in the presence of external charges is 
the correlator 

The matrix element in (^) has the standard phase space path integral 
expression |54], so that 



e 



(100) 

where the time interval is r G [0, 1/T], the spatial integral in the action is 
taken over S'^, the canonical momentum vr has open boundary conditions and 
the gauge field has the boundary condition which is periodic up to a twist by 
a gauge transformation, 

a(l/T, x) = a(0, x) - WAq{x)/T (101) 
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We have denoted the fluctuating part of the magnetic field as 6 = V x 
a. Equation ( |10U| ) gives the effective action up to an overall temperature 
dependent but Aq independent constant. 

The canonical momentum can be integrated in order to present the func- 
tional integral in configuration space, up to an irrelevant temperature depen- 
dent factor, 

(102) 

It is convenient to untwist the boundary condition using the change of 
variables, 

a(r, x) a{r, x) — tVAq^x) (103) 

with the result that 

Bm 

■ J [da] e~ ''^ '^^ (i<j'+ii?i,+i62+i|sxS-«Aob) ^^^^^ 

The Gaussian integral over a, can now easily be done. If we choose Bm to be 
a constant magnetic field on the sphere, the integral over a{x) in (|104|) yields 



2 



e 



1 



2T(47ri?2) {47tR'^)T Js^ 



E \Aoil,m)\\l{l + l)/R' + K')\ (105) 

where we have dropped an irrelevant infinite constant, recalled the fact that 
Js2 = g = AttR'^Bm and we have expanded Ao{x) in spherical harmonics, 

Ao{x)=J2Ml,m)Yimix) (106) 

Im 

I > are integers of the usual angular momentum spectrum, the spectrum 
of the laplacian — being /(/ + 1)-R^, and m = —I, + for each /. 

Notice that the non-zero modes of the gauge function x{^) are governed by 
a massive euclidean free field theory. The zero modes, on the other hand are 
coupled to the monopole moments which must be summed. 

If we recall the monopole quantization condition, g = 27m/ e, the summa- 
tion over monopoles is 

where 

Using the Poisson resummation formula, the sum can be presented as 

/(47ri?2)e2T^ / {A^R^)e'T f k \2\ 
^^^Eexp(-^^(^Ao + n) j (108) 
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which exphcitly exhibits periodicity in Aq, 



eT 

Ao^Ao + — (109) 

K 

The summation over monopoles effectively represents this symmetry in the 
Villain form 



g-Scff [Ao] 



This action is identical (although one dimension higher) to the action that 



was found for the Schwinger model in ref. ||T9|. In that case, we argued that 
the Z symmetry is always broken in the infinite volume limit. This symmetry 
breaking is interpreted as screening, similar to that which occurs in a Higgs 
phase, rather than deconfinement. We conjecture that a topologically massive 
gauge theory screens, rather than deconfines. We shall derive support for 
this conjecture from the variational calculation in Section IV where we find 
indications of a rapid change of the behavior of the system between what 
we here interpret as a screening phase and what we would properly call a 
deconfined phase. 

In the next subsection we shall examine the consequences of this symmetry 
of the effective action ( |11U| ) and we will show that it is also spontaneously 
broken in the infinite volume limit. 

3.3 Spontaneous breaking of Z symmetry 

The Polyakov loop operator transforms under ( p.U9| ) as 

^ineAo(x)/T ^ ^ine^ / k ^ineAo{x)/T _ ^2nin-p/q ^ineAo{x)/T (HI) 

where we have defined (as in ( |5^ ) 

2 

e p 

This symmetry has implications for the correlator of Polyakov loop operators, 

Z[T, (e„ Xi)] = (jl e'^^^'^^-^y^^ (112) 

which depend on the ratio p/q. 

• p/q is irrational. Then the correlator ( p.l2| ) vanishes unless ^iCi = 0. 

• p/q is rational. Then the correlator (|112|) vanishes unless J2j = 
e ■ integer ■ q. Since the consistency with the monopole bundle requires 
that the charges are quantized in units of e, so J2j = integer ■ e, the 
Z symmetry here is actually a subgroup of Z, Zq, the additive group of 
the integers modulo q. 

• If p is an integer and q = 1 the correlator ( |112| ) is unrestricted. 
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It is interesting to observe that only charges which are integer multiples of 
qe are allowed on the sphere. This is no surprise as it is a direct consequence of 
Gauss' law. The integral of Gauss' law together with the Dirac quantization 
condition yields the constraint charge = qe ■ integer . Thus, only charges 
which are integer multiples of q are consistent with the Gauss' law constraint 
on the sphere. The Zq symmetry enforces this "topological constraint" . 

The Zq-symmetry is invariably broken in the infinite volume limit. To see 
this, we consider the two-point correlator of loop operators 

^ieAo{x)/T ^-ieAo{0)/T^ ^ . (^-^{x\ ^ ^ |0) j (113) 

If we take the limit as the volume goes to infinity, followed by the limit as 
the separation of the points in the correlator goes to infinity, the two-point 
function approaches a non-zero constant, 

lim lim ( e^"^°(^)/^ e-^"^°(°)/^\ = const. (114) 

|sH0Oii2-»0O \ / 

This implies that the Zq symmetry is spontaneously broken in the infinite 
volume limit. 

This leaves us with the correct conclusion that the topologically massive 
gauge theory is not a confining theory. In fact its electrostatic interactions are 
short-ranged and Yukawa-like. Their large distance fall-off is governed by the 
inverse topological mass. 

Thus, for all practical purposes, the topological mass of the photon 
contributes a mass term to the effective action for Aq. The fact that this 
mass term is periodic when the volume is finite is irrelevant in the infinite 
volume limit. If we couple topologically massive QED to matter fields, we 
would expect that, in the limit where the matter field masses are large, the 
effective field theory is the massive sine-Gordon theory, 

ST^^iAo] = I d^x I^V^o • V^o + |;^o - ^ e-/^cos(e^o/r)} 

(115) 

In the next Section, we shall study this model using a variational approach. 



4 Variational approach 

In this section, we shall discuss a variational approach to the problem of 
showing the existence of a phase transition in the one-loop effective theory. 
Wc have argued that the effective field theory where the phase transition can 
be studied is the sine-Gordon theory with an additional mass term for the 
boson, 

^efr = /ci'x|^V0-V</)+^02-^cos/?</)| (116) 

Here, the mass term for the boson is the topological photon mass, /i — k. 
Also, in the limit where the fermion masses of QED are much greater than the 
temperature and charge squared, we have /3 = e/y/T and a — e^m e~"^/-^/7r. 
The ultraviolet cutoff is the fermion mass, m. 
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As we have discussed in the previous sections, the mass term for the boson 
should really be a periodic one. However, quantum effects invariably break 
the translation symmetry for such a Bose field in the limit where the volume 
is infinite. We shall therefore ignore the translation symmetry. 

In Ref. the behavior of the model ( |116| ) was investigated around the 



point = 47r, and a phase transition at the corresponding critical temperature 
Tc = e^/47r was suggested. In subsequent work on the massless model it 
was however found that the divergencies at (3"^ = Air can be handled within the 
renormalization scheme, and the behavior is smooth. Accordingly, we shall 
concentrate on studying the model only in the neighborhood of = Svr. 

In Section II we reviewed the argument of ref. |]T7| that, when the topo- 
logical mass is zero, there is a Berezinskii-Kosterlitz-Thouless confinement- 
deconfinement transition at the temperature Tc[e^, m] = e^/87r(l + ...C(e^/m)), 
corresponding to ~ Sir for the sine-Gordon parameter. When there is 
a topological mass, the system always fails to attain its asymptotic scaling 
regime, and the flow diagram cannot be that of a BKT transition. In fact, at 
a scale of the order of the inverse topological mass there is a crossover from a 
behavior governed by the long ranged Coulomb interaction to one governed by 
a finite range Yukawa interaction. There is no confinement in the strict sense, 
since the Coulomb interaction cuts off at a given distance. However, when 
the topological mass is small, much smaller than any of the other dimensional 
parameters (particularly the confining scale which is governed by e^), for low 
temperature and spatial scales much less than the inverse topological mass 
and of the order of the confinement scale, the physical behavior of the sys- 
tem should be much like a confining one with a confining interaction between 
oppositely charged particles. One expects that, at or near the confinement- 
deconfinement temperature, a drastic change in the properties of the system 
takes place, though not a true phase transition. Such a rapid change is be- 
tween a quasi- confining and a de-confined behavior with the electric mass of 
the photon moving from the value of the topological mass, in the low temper- 
ature quasi-confined phase, to the value of the much larger Debye screening 
mass, in the de-confined phase. 

If one imagines separating a test particle-antiparticle pair at short distances 
the potential energy should increase with distance just as it does in a strictly 
confining system. If this increasing energy gets large enough it can produce 
a pair of dynamical charges which partially screen the charges of the source. 
Thus, at short distances, separating a particle-antiparticle pair should produce 
mesons, rather than free charged particles. However, at distances larger than 
the inverse of the topological mass, charged particles should behave as free 
particles. We shall interpret this as a screening (as distinct from deconfined) 
phase analogous to what happens in the Schwinger model or massless two- 
dimensional QCD. 

If we then increase the topological mass to the confining scale, the system 
should go continuously to one which deconfines at all scales. We shall 
qualitatively study these behaviours by means of a variational approach in 
the next subsections. 
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4.1 Jensen's inequality 

In statistical mechanics, a variational approach uses Jensen's inequality. First, 
we shall give a brief review of this inequality and its derivation. Consider a 
statistical system with Hamiltonian Hi which we want to study the statistical 
mechanics of, but are unable to solve for the sum over states to obtain 
the partition function or the correlators exactly. We consider another test 
Hamiltonian, Hq which contains some parameters and with which we can solve 
for the partition function and correlation functions of observables analytically. 
Then, consider the Hamiltonian which interpolates linearly between them. 



H^^XHi + {l-X)Ho (117) 
W{X) = -In^ e-^^(^) (118) 



^^^^^ - (H H) 
Hi-Hof)+{Hi-Ho)l < (119) 



and the free energy 
It has the properties 

~dXP~ " /A 

where {■■■)x is the expectation value in the ensemble with Hamiltonian H,. 
Since the curvature of 1^(A) as a function of A is always less than or equal to 
zero, 1^(A) obeys the inequahty 

W{X)<W{0) + X-(^^^'^^ ^ (120) 

which, evaluated at A = 1 is Jensen's inequality 

W(l)<W(0) + {Hi-Ho)o (121) 

This establishes an upper bound on the free energy of the system of interest 
by the system in the variational ansatz which can be optimized by adjusting 
the parameters of the variational ansatz. The bound is saturated only when 
the ensembles are identical (Hi — Hq). 

4.2 Variational computations 

We are faced with the problem of computing the thermodynamics of the system 
with Hamiltonian function 

Hi = J d'x I \d,<pd,<P + y 0' - ^ cos (122) 

We shall study this system variationally by beginning with the test ensemble 
governed by the quadratic Hamiltonian 

Ho^J d'xd'y^mai^ - yf)<t>{y) (123) 
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with an ultraviolet cutoff A, which is of order the mass of the matter field in the 
original model. The Gaussian functional integrals appearing on the right-hand 
side of Jensen's inequality in this case can easily be done with result 

V<ir|.nf^V4^#-l|-^e-f/V- . ,124) 



V - 2J { \ J J /32 

Here the measure for the cutoff momentum integration is 



e{A'-p') (125) 



and has been extracted a factor of the spatial volume to get the free energy 
density. 

In order to optimize the ansatz, we take the variational derivative of the 
right hand side to obtain 



= ap'r'i^ip') -p^-fi^-ae'^l"^^-"] (126) 



The equation for the minimum is solved by 

e(/)=p' + M2 (127) 
where the mass parameter satisfies the equation 

We substitute (|127| ) into ( |124|) and integrate to obtain 



W< — \A^\n{ +/i2ln — — -Svr 
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V - 8tt \ V A2 y ^ V / f3^ \A^ + My 

(129) 

The regime that we are interested in is where the cutoff, (the electron mass) is 
very large and the other dimensional parameters n and a are very small. We 
look for a minimum of the potential ( p.29| ) where M is of the same order of 
magnitude as fi and a. This corresponds to seeking a solution of sine-Gordon 
theory which is consistently renormalized as a relativistic quantum field theory. 

4.2.1 Sine-Gordon theory (yU^ = q, A ^ oo) 



It is interesting to explore the minima of (|129|) in the case of pure sine-Gordon 
theory, when /i^ = 0. This can be done using a slightly different variational 
method by Coleman Q . 

We must seek a solution of the equation for the variational mass where 
<< A^. To do this, we must first renormalize the bare coefficient of the 
cosine term in the action, 

a = aRi — \ (130) 
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where a is an arbitrary mass scale which accompanies renormahzation. Then, 
up to an infinite, M independent constant, the variational free energy is 

Clearly, when /J^/Svr < 1, the point M = is a local maximum of the 
potential. There is a minimum at M 7^ 0, determined by the finite dimensional 
parameters a and a/j and by f3. On the other hand, when > 1, there is 

a local minimum of the potential at M = but for large M, it is unbounded 
from below. This means that M = is an unstable state and the value of M 
at the true minimum of the potential is of order the (infinite) cutoff. 

Thus, we regain Coleman's result [Q. Without the renormahzation of the 



parameter, which was shown to be necessary in ref. [^, the pure sine- 
Gordon model has a phase transition at the point = Stt from a phase where 
the theory is approximately solved by a boson with finite mass to one where 
the boson has infinite mass, of order of the cutoff. 



4.2.2 Sine-Gordon theory: parity invariant QED (fi^ = 0, A finite ) 

Let us now analyze the case in which the cutoff is finite and of the order of 
the electron mass m. We can then analyze directly the potential ( |129D with 
A ~ m and /i = 0. It can be easily seen that, independently on the value 
of A, the potential ( p.29| ) for /^^/Svr < 1 behaves as before, i.e. at the point 
M = has a local maximum and an absolute minimum at M 7^ 0. For 
P'^/SiT > 1, however the minimum at M = is now an absolute minimum for 
any M < A, the potential is always bounded below. M = describes a stable 
state, so that the transition is between a phase characterized by the behavior 
of a two-dimensional massless boson to one charcterized by a massive boson. 

This can be interpreted as a change in the symmetry, since a massless boson 
has the field translation symmetry, 0(x) — > 0(a;) + const, whereas a massive 
boson does not. If we translate the parameters of the sine-Gordon theory 
into those of the effective action for QED (as in the discussion after equation 
(|116|) ), the phase transition occurs at the critical temperature Tc = e^/Svr. 
Above this temperature, the boson has a mass, corresponding to deconfined 
phase which screens electric charges by virtue of the "Debye" mass M (region 
III in Fig.l), whereas below this temperature the boson is massless (region 
I in Fig.l), corresponding to a confining phase which cannot screen the long 
ranged electric fields of incommensurate charges. 

4.2.3 Sine-Gordon theory with a mass (/x^ 7^ 0, A ^ 00) 

Let us now consider the case of massive sine-Gordon theory with /i^ > 0. 
The variational potential is made finite by the same renormahzation of the 
parameter a as in the massless case. The variational free energy is then 

l«,._L(M^_,.,n(M>V8.^(^)'''") (132, 

If < Stt, the minimum of the potential occurs at a finite value of M^. This 
value depends crucially also on the topological mass fi and for fi very small 
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occurs approximately at the same value of the = case, for sufficiently large 
it occurs at fi. As in massless sine-Gordon theory, without a renormalization 
of /?^, the potential is unbounded below if > Svr, independently on the value 
of n- This implies that the global minimum of the potential is of order the 
ultraviolet cutoff. 



4.2.4 Sine-Gordon theory with a mass: topologically massive QED 
(/i^ ^ 0, A finite) 

For /i very small and (3"^ < Sir (high temperature) the potential ( p.29|) has an 
absolute minimum at a finite value of M much larger than n, independently 
on the value of the cutoff A . For > Svr (low temperature) the absolute 
minimum rapidly moves to the small value M ~ /i, so that the transition is 
between two distinct massive behaviors. For a very tiny fi the crossover in 
the BKT should arise only at very large scales. The behavior of the system 
in the two phases should be quite different. In the low temperature region all 
the components of the gauge field have a small mass (the topological mass), 
and this phase should be very much like the Higgs phase of the Schwinger 
model |]l9l. In the high temperature region, the electric mass grows to the 



much bigger Debye mass. This region is essentially similar to the plasma 
phase of the parity invariant theory. Moreover, the Debye screening mass 
is temperature dependent, whereas the topological mass does not depend on 
temperature. 

It is interesting that this transition occurs even when there is an explicit 
mass term in the action, provided that the mass term is sufficently small. 
Although this transition is in a sense associated with vortex binding-unbinding, 
just as it is in the sine-Gordon theory which describes the Coulomb gas, it is 
not a confinement-deconfinement transition in the strict sense, since the Z 
symmetry is broken in both phases. 

When the topological mass fi is large, the minimum of the potential is 
always at /i regardless of the value of so that the crossover arises at short 
scales (~ l//i) and the fiow diagram is completely destroyed. No drastic change 
in the behavior of the system should be observed as the temperature is changed 
in the neighborough of Tc = c^/Stt. 



5 Discussion 

In this paper, we have shown that deconfinement in finite temperature QED 
can be characterized as breaking of a certain global discrete symmetry. We 
have also shown that a confinement-deconfinement transition takes place in 
parity invariant 2+1-dimensional QED, at least in the regime where the 
electron mass is large. 

In a sense, the latter fact is no surprise. When the electromagnetic coupling 
is small compared to the electron mass m so that vacuum fiuctuations 
are suppressed, and when the temperature is also smaller than the electron 
mass, the thermal state is to a good approximation a dilute two-dimensional 
neutral Coulomb gas of thermally excited electrons and positrons. It is well 
known that this Coulomb gas, even when very dilute, has a BKT transition. 
Below the transition temperature, the electrons and positrons are bound 
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into pairs. Above the transition temperature, electrons and positrons are 
approximately free particles. One physical prediction which can be deduced 
from the presence of the BKT transition is the universal property of the phase 
transition associated with the bulk modulus of spin waves in the gapless, 
Kosterlitz-Thouless phase. The exponents in the correlators of Polyakov loop 
operators in that phase are determined by a single correlation length divided 
by the temperature. 

There is a similar picture of the topologically massive theory. If the 
topological mass is small, at weak coupling and temperatures somewhat less 
than the electron mass, the thermal state is to a good approximation that of 
the dilute Yukawa gas. There is a rapid change between a quasi-confining and 
a deconfined behavior close to the BKT critical point of the parity invariant 
theory. 

In both cases, our analysis is only valid where the electron mass is large. 
We expect that the phase transition, or more correctly the BKT line of phase 
transitions, persists for some time as we lower the electron mass or raise the 
electric charge. However, the resulting strong coupling regime is out of the 
domain of validity of our analysis. 

An interesting, experimentally testable consequence of the breaking of the 
Z symmetry is the existence of domain walls. For the existence of such, it 
is enough to have field configurations that are non-contractable loops in the 
f/(l) group. The question of domain walls has been addressed in gluodynamics 
p3| , |3^ and massless QED [|1^. Whether they correspond to real Minkowski 
space objects, is an interesting open question. It is instructive to note that 
upon refermionizing our effective sine-Gordon model (|55D along the lines of 



44| , domain walls in sine-Gordon theory correspond to worldlines of fermions 
in the resulting Thirring model. 

Phase transitions analogous to the one which we conjecture to exist in 
topologically massive QED have been studied experimentally in quasi two- 
dimensional condensed matter systems, particularly charged vortex arrays in 
superconducting films and have also been useful in theoretical work on 
high Tc superconductivity where a type of "smoothed" BKT transition 
is discussed. The experimental study of analog systems in condensed matter 
physics could help to resolve some of the theoretical questions raised by our 
current work. 



Acknowledgement 

We thank I. Affleck, A. D'Adda, R. Ferrari, A. Kovner, I. Kogan, A. Polyakov, 
P. C. E. Stamp, L.C.R. Wijewardhana and N. Weiss for helpful discussions. 
G. Semenoff and O. Tirkkonen acknowledge the kind hospitality of the Dipar- 
timento di Fisica, Universita di Perugia and Universita di Trento where part 
of this work was done. 



References 

[1] W. Marciano and H. Pagels, Phys. Rep. 36, 137 (1978). 



28 



[2] For a modern review, see A. Di Giacomo, Lecutes at CXXX Course 
of the International School of Physics E. Fermi on "Selected topics in 
nonperturbative QCD", Varenna, June 1995., hepth/9603029. 

[3] For reviews, discussion and further references, see J. Swieca in "Field 
Theoretical Methods in Particle Physics", W. Ruhl, Editor, Plenum, New 
york, 1979, pg. 325 and E. Abdalla, M. Abdalla and K. Rothe, "Non- 
perturbarive methods in two-dimensional quantum field theory" (World 
Scientific, Singapore, 1991). 

[4] D. Gross, I. Klebanov, A. Matytsin and A. Smilga, Princeton preprint 
|hep-th/95TTT0^ . 

[5] A.M. Polyakov, Phys. Lett. 72B, 477 (1978). 

[6] L. Susskind, Phys. Rev. D20, 2610 (1979). 

[7] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 

[8] L. McLerran and B. Svetitsky, Phys. Lett. 98B, 195 (1981); Phys. Rev. 
D24, 450 (1981). 

[9] J. Kutti, J. Polonyi and K. Szlachanyi, Phys. Lett. 98B, 199 (1981). 
[10] J. Engels, F. Karsch, H. Satz and I. Montvay, Phys. Lett. lOlB (1981). 
[11] B. Svetitsky and L. G. Yaffe, Nucl. Phys. B210, 423 (1982). 
[12] N. Weiss, Phys. Rev. D44, 75 (1981). 

[13] M. Billo, M. Caselle, A. D'Adda, L. Magnea, S. Panzeri, Nucl. Phys. 
B435, 172 (1995). 

[14] M. Billo, M. Caselle, A. D'Adda, S. Panzeri, Torino preprint DFIT-69-95, 
|hep-lat/960 10201 . 

[15] N. Weiss, Phys. Rev. D35, 2495 (1987); A. Roberge and N. Weiss, Nucl. 
Phys. B275, 734 (1986). 

[16] A. V. Smilga, Ann. Phys. (N.Y.) 234, 1 (1994). 

[17] G. Grignani, G. Semenoff and P. Sodano, Perugia and UBC preprint 
(1995), |hep-th/9504105| . 

[18] E. Braaten and A. Nieto, Phys. Rev. Lett. 74, 3530 (1995). 

[19] G. Grignani, G. Semenoff, P. Sodano and O. Tirkkonen, Perugia and UBC 
preprint (1995), |hep-th/951lTT0| . 

[20] K. Johnson, M. Baker and R. Willey, Phys. Rev. Lett. 11, 518 (1963); 
Phys. Rev. 136 (1964), 1964; V. Miransky, Nuov. Cim. 90A, 149 (1985). 

[21] H. Hansson, H.B. Nielsen and I. Zahed, Nucl. Phys. B451, 162 (1995). 

[22] I. Kogan, Phys. Rev. D49, 6799 (1984). 



29 



[23] D. Kutasov, Nucl. Phys. B414, 33 (1994). 

[24] G. Bhanot, K. Demeterfi and 1. Klebanov, Phys. Rev. D48, 4980 (1993); 
K. Demeterfi, G. Bhanot and I. Klebanov, Nucl. Phys. B418, 15 (1994). 

[25] I. Kogan and A. Zhitnitsky, UBC and Oxford preprint (1995), |hep- 
I ph/950932'2| . 

[26] F. Lenz, M. Shifman and M. Thies, Phys. Rev. D51, 7060 (1995). 
[27] D. Sen, Phys. Rev. D41, 1227 (1990). 

[28] T. Appelquist, D. Nash and L. C. R. Wijewardhana, Phys. Rev. Lett. 60, 
2575 (1988). 

[29] V. Berezinksii, Sov. Phys. JETP 34, 610 (1972) [Zh. Eksp. Teor. Fiz. 61, 
1144 (1972)]. 

[30] J. M. Kosterlitz and D. J. Thouless, J. Phys C: Sohd St. Phys. 6, 1181 
(1973). 

[31] R. Jackiw and S. Templeton, Phys. Rev. D23, 2291 (1981). 

[32] W. Siegel, Nucl. Phys. B156, 135 (1979); J. Schonfeld, Nucl. Phys. B185, 
157 (1981); S. Deser, R. Jackiw and S. Templeton, Phys. Rev. Lett. 48, 975 
(1982); Ann. Phys. 140, 372 (1982). 

[33] T. Bhattacharya, A. Gocksch, C. Korthals-Altes and R. D.Pisarski, Phys. 
Rev. Lett. 66, 998 (1991). 

[34] V. Dixit and M. G. Ogilvie, Phys. Lett. B269, 353 (1991); J. Ignatius, K. 
Kajantie and K. Rummukainen, Phys. Rev. Lett. 68, 737 (1992). 

[35] J. Kiskis, Phys. Rev. D51, 3781, (1995); preprint |hep-lat/951002£ . 

[36] A. Niemi and G. Semenoff, Phys. Rev. Lett. 51, 2077 (1983). 

[37] N. Redlich, Phys. Rev. Lett. 52, 18 (1983); Phys. Rev. D29, 2366 (1984). 

[38] T. Appelquist and R. Pisarski, Phys. Rev. D23, 2305 (1981). 

[39] T. M. Rice, Phys. Rev. 140, 1889 (1966); J. W. Kane and L. P. Kadanoff, 
Phys. Rev. 155, 80 (1967); F. J. Wegner, Z. Phys. 206, 465 (1967); J. Zittatrz, 
Z. Phys. B23, 55 (1976). 

[40] For an up to date discussion, see N. Goldenfeld, "Lectures on Phase 
Transitions and the Renormalization Group" , (Addison- Wesley, New York, 
1992). 

[41] R. Bahan, C. Itzykson, G. Parisi and J. B. Zuber, Phys. Rev. D17, 1041 
(1978). 

[42] I. S. Gradshteyn and L M. Ryzhik, "Table of Integrals Series and Prod- 
ucts" (Academic Press, Inc., London, 1980). 



30 



[43] L. Lewin, "Dilogarithms and Associated Functions" (Macdonald, London, 
1958). 

[44] S. Coleman, Phys. Rev. Dll, 2088 (1975). 

[45] S. T. Chui and P. A. Lee, Phys. Rev. Lett 35, 315 (1975). 

[46] S. Samuel, Phys. Rev. D18, 1916 (1978) 

[47] P.B. Wiegmann, J. Phys. C7, 1583 (1978). 

[48] J. V. Jose, L. P. Kadanoff, S. Kirkpatrick and D. R. Nelson, Phys. Rev. 
B16, 1217 (1978). 

[49] D. J. Amit, Y. Y. Goldschmidt and G. Grinstein, J. Phys. A13, 585 (1980); 
D. J. Amit, S. Elitzur and E. Rabinovici and R. Savit, Nucl. Phys. B210, 
69 (1982). 

[50] W. Chen, M. Dobroliubov, G. Semenoff, Phys. Rev. D46, 1223 (1992). 
[51] V. Belyaev, I. Kogan, G. Semenoff, N. Weiss, Phys. Lett. B277, 331 (1992). 
[52] M. Bergeron and G. Semenoff, UBC preprint (1993), |hep-th/ 93060^ . 
[53] O. Alvarez, Comm. Math. Phys. 100, 279 (1985). 

[54] L. Schulman, "Techniques and application of path integration" (Wiley, 
New York, 1981). 

[55] J. Villain, J. Physique 36, 581 (1975). 

[56] J. Frdhhch, Commun. Math. Phys., 47, 233 (1976). 

[57] A. F. Hebard and A. T. Fiory, Physica B 109, 1637 (1982); A. M. Kadin, 
K. Epstein and A. M. Goldman, Phys. Rev. B27, 6691 (1983); A. T. Fiory, 
A. F. Hebard and W. L Glaberson, Phys. Rev. B28, 5075 (1983); J. C. 
Garland and H. J. Lee, Phys. Rev. B36, 3638 (1987). 

[58] P. C. E. Stamp, Phys. Rev. B38, 2847 (1988). 



31 




87t e/T 



Figure 1. 

Renormalization flow diagram for tlie BKT transition. Tlie arrows denote flow 
toward the ultraviolet. Region I is the confining phase whereas regions II and 
III are deconfined. Region III is asymptotically free whereas in region I there 
is a line of infrared stable fixed points which represent c = 1 conformal field 
theories. The separatrix between regions 1 and 11 is the line of BKT phase 
transitions. The critical behaviour of the system at the latter phase transition 
is that of an SU(2) Wess-Zumino-Witten model. 
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